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Second order QCD corrections to jet production at hadron colliders: 

the all-gluon contribution 
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We report the calculation of next-to-next-to-leading order (NNLO) QCD corrections in the purely 
gluonic channel to dijet production and related observables at hadron colliders. Our result represents 
the first NNLO calculation of a massless jet observable at hadron colliders, and opens the path 
towards precision QCD phenomenology with the LHC. 

PACS numbers: 13.87.Ce,12.38Bx 



Single inclusive jet and dijet observables are the most 
fundamental QCD processes measured at hadron collid- 
ers. They probe the basic parton-parton scattering in 
2 — > 2 kinematics, and thus allow for a determination of 
the parton distribution functions in the proton and for 
a direct probe of the strong coupling constant a s up to 
the highest energy scales that can be attained in collider 
experiments. 

Precision measurements of single jet and dijet cross 
sections have been performed by CDF [l[ and DO Q at 
the Tevatron and by ATLAS Q and CMS 0| at the LHC. 
The Tevatron data are included in nearly all global fits of 
parton distributions, where they provide crucial informa- 
tion on the gluon content of the proton, and have been 
used to determine the strong coupling constant 0. 

Theoretical predictions for these observables are accu- 
rate to ncxt-to-leading order (NLO) in QCD ftjl-fic]| and 
the electroweak theory The estimated uncertainty 

from missing higher order corrections on the NLO QCD 
predictions is substantially larger than the experimen- 
tal errors on single jet and dijet data, and is thus the 
dominant source of error in the determination of a s . A 
consistent inclusion of jet data in global fits of parton dis- 
tributions is also feasible only to NLO. These theoretical 
limitations to precision phenomenology provide a very 
strong motivation for computing next-to-next-to-leading 
order (NNLO) corrections to jet production at hadron 
colliders. 

At this perturbative order, three types of parton-level 
processes contribute to jet production: the two-loop vir- 
tual corrections to the basic 2 — » 2 process [l2j, the 
one-loop virtual corrections to the single real radiation 
2 — >• 3 process 13j and the double real radiation 2 — s- 4 
process at tree-level [3]. Each contribution is infrared 
divergent, and only their sum yields a finite and mean- 
ingful result. After ultraviolet renormalization, both vir- 
tual contributions contain explicit infrared singularities, 
which are compensated by infrared singularities from sin- 
gle or double real radiation. These become explicit only 
after integrating out the real radiation contributions over 
the phase space relevant to single jet or dijet production. 



This interplay with the jet definition complicates the ex- 
traction of infrared singularities from the real radiation 
process. It is typically done by subtracting an infrared 
approximation from the corresponding matrix elements. 
These infrared subtraction terms are sufficiently simple 
to be integrated analytically, such that they can be com- 
bined with the virtual contributions to obtain the cancel- 
lation of all infrared singularities. Several generic meth- 
ods for the construction of subtraction terms are available 
at NLO EMll- 



21] 



The development of subtraction methods for NNLO 
calculations is a very active field of research. Up to 
now, various methods were constructed and applied to 
specific NNLO calculations of exclusive observables: sec- 
tor decomposition [3] applied to Higgs production 19 1 
and vector boson production [20| : g^-subtraction 
to Higgs production [22j] . vector boson production 
associated F-ff-production [24|, photon pair produc- 
tion 12511 and top quark decay 26] j antenna subtrac- 
tion [23 to three-jet production [23, |29j and related event 
shapes |3fll3l"l] in e + e~ annihilation; and sector-improved 
residue subtraction 32] to top quark pair production [33| . 

The antenna subtraction method 27, 34, 35| constructs 
subtraction terms from antenna functions which encap- 
sulate all unresolved radiation in between a pair of hard 
radiator partons. At NNLO, antenna functions with up 
to two unresolved partons at tree level and one unresolved 
parton at one loop are required. For hadron collider ob- 



servables, one [36| or both [37|, l38j radiator partons can 



be in the initial state. 

For the NNLO all-gluon contribution to jet produc- 
tion at hadron colliders, the antenna subtraction terms 
were constructed for the tree-level double real radiation 
process in |39j and for the one-loop single real radiation 
process in [40j. These subtraction terms were integrated 
and combined 41[ with the relevant two-loop matrix el- 
ements and parton distributions, resulting in a full can- 
cellation of infrared poles. We have now implemented 
these terms into a parton-level event generator, which 
can compute the all-gluon contribution to any infrared- 
safe observable related to dijet final states at hadron col- 



2 



liders to NNLO accuracy. The program consists of three 
integration channels: 



d<7 



gg.NNLO 



d* 4 



d* 3 



d* 2 



[da 



RR da S 



gg,NNLO 



gg,NNLO\ 



[ d °fg,NNLO - d& 



l aa gg,NNLO 



gg,NNLo] 



where each of the square brackets is finite and well be- 
haved in the infrared singular regions. For the all-gluons 
channel, the construction of the three subtraction terms 
ij^NNLO was described in Refs. 



39H41|. 



do-; 

In the three-parton and four-parton channel, the phase 
space has been decomposed into multiple wedges (6 
three-parton wedges and 30 four-parton wedges), each 
containing only a subset of possible infrared singular con- 
tributions. Inside each wedge, the generation of multiple 
phase space configurations related by angular rotation of 
unresolved pairs of particles around their common mo- 
mentum axis ensures a local convergence of the antenna 
subtraction term to the relevant matrix element. Owing 
to the symmetry properties of the all-gluon final state, 
many wedges yield identical contributions, thereby al- 
lowing a substantial speed-up of their evaluation. 

Jets in hadronic collisions can be produced through 
a variety of different partonic subprocesses, and the all- 
gluon process is only one of them. Our results on this 
process can therefore not be directly compared with ex- 
perimental data. The all-gluon process does however al- 
low to establish the calculational method, and to qualify 
the potential impact of NNLO corrections on jet observ- 
ables. It should be noted that the NLO corrections to 
hadronic two- and three-jet production were also first 
derived in the all-gluon channel 142 -|44 j , well before full 
results could be completed [10, El- In both cases, the 
all-gluon results were extremely vital both for establish- 
ing the methodology and for assessing the infrared sensi- 
tivity of different jet algorithms [44| . 

Our numerical studies for proton-proton collisions at 
centre-of-mass energy \fs = 8 TeV concern the single 
jet inclusive cross section (where every identified jet in 
an event that passes the selection cuts contributes, such 
that a single event potentially enters the distributions 
multiple times) and the two-jet exclusive cross section 
(where events with exactly two identified jets contribute). 

Jets are identified using the anti-fcr algorithm with res- 
olution parameter R = 0.7. Jets are accepted at central 
rapidity \y\ < 4.4, and ordered in transverse momentum. 
An event is retained if the leading jet has pn > 80 GeV. 
For the dijet invariant mass distribution, a second jet 
must be observed with pti > 60 GeV. 

All calculations are carried out with the 
MSTW08NNLO gluon distribution function £3, 
including the evaluation of the LO and NLO contri- 
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FIG. 1: Inclusive jet transverse energy distribution, da/dpT, 
for jets constructed with the anti-fcr algorithm with R — 0.7 
and with p T > 80 GeV, \y\ < 4.4 and y/s = 8 TeV at NNLO 
(blue), NLO (red) and LO (dark-green). The lower panel 
shows the ratios of NNLO, NLO and LO cross sections. 



butions [47j. This choice of parameters allows us to 
quantify the size of the genuine NNLO contributions 
to the parton-level subprocess. Factorization and 
rcnormalization scales {jip and fj-u) are chosen dynami- 
cally on an event-by-event basis. As default value, we 
set fip = /i_R = /i an d se t A* equal to the transverse 
momentum of the leading jet so that ^ — pxi- 
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FIG. 2: Scale dependence of the inclusive jet cross section for 
pp collisions at = 8 TeV for the anti-fcr algorithm with 
R = 0.7 and with \y\ < 4.4 and 80 GeV < p T < 97 GeV at 
NNLO (blue), NLO (red) and LO (green). 

In Fig. Q] we present the inclusive jet cross section for 
the anti-/cT algorithm with R — 0.7 and with px > 
80 GeV, \y\ < 4.4 as a function of the jet px at LO, 
NLO and NNLO, for the central scale choice /i = pn- 
The NNLO/NLO k- factor shows the size of the higher 
order NNLO effect to the cross section in each bin with 



3 



respect to the NLO calculation. For this scale choice we 
see that the NNLO/NLO A:- factor is approximately flat 
across the px range corresponding to a 15-25% increase 
compared to the NLO cross section. 

One of the main motivations for computing the NNLO 
QCD corrections is to reduce the scale uncertainty in the 
theoretical prediction. This is illustrated in Fig. [2] for the 
single jet inclusive cross section for jets with \y\ < 4.4 and 
80 GeV < p T < 97 GeV. We see that the scale depen- 
dence of the cross section at NNLO is vastly reduced. 
The scale dependence of other px and y slices is also 
reduced. 

To illustrate the range of observables that can be stud- 
ied with our computation we show in Fig. |3] the inclusive 
jet cross section in double-differential form in jet pr and 
rapidity bins at NNLO. The pt range is divided into 16 
jet-pr bins and seven rapidity intervals over the range 
0.0-4.4 covering central and forward jets. 
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FIG. 4: Double differential fc-factors for pr > 80 GeV and 
three \y\ slices: \y\ < 0.3, 0.3 < \y\ < 0.8 and 0.8 < \y\ < 1.2. 
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FIG. 3: The doubly differential inclusive jet transverse energy 
distribution, d 2 a/dpTd\y\, at y^s = 8 TeV for the anti-fcr 
algorithm with R = 0.7 and for Et > 80 GeV and various \y\ 
slices. 

Fig. 0] shows the double-differential fc-factors for the 
distribution in Fig. [3] for three rapidity slices: \y\ < 0.3, 
0.3 < \y\ < 0.8 and 0.8 < \y\ < 1.2. We observe that 
the NNLO correction increases the cross section between 



25% at low px to 12% at high px with respect to the 
NLO calculation and this behaviour is similar for all three 
rapidity slices. 
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FIG. 5: Exclusive dijet invariant mass distribution, 
da/dmjjdy", at yfs = 8 TeV for y* < 0.5 with pri > 80 GeV, 
PT2 > 60 GeV and \ yi \, \y 2 \ < 4.4 at NNLO (blue), NLO 
(red) and LO (dark-green). The lower panel shows the ratios 
of NNLO, NLO and LO cross sections. 

As a final observable, we computed the dijet cross 
section as a function of the dijet mass at NNLO. This 
is shown in Fig. [5] for the scale choice /i = 2p^i to- 
gether with the LO and NLO results. The dijet mass 
is computed from the two jets with the highest pt and 
1 2/1 1 , I2/2I < 4.4 with y*, defined as half the rapidity dif- 
ference of the two leading jets y* = \yi — x/2 1 / 2 < 0.5. We 
see that the NNLO/NLO fc-factor is approximately flat 
across the rtijj range corresponding to a 15-20% increase 
compared to the NLO cross section. 

In conclusion, we have described the first calculation of 
the fully differential inclusive jet and dijet cross sections 
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at hadron colliders at NNLO in the strong coupling con- 
stant using the new parton-level generator NNLOJET. 
We have considered the NNLO QCD corrections from 
the purely gluonic channel at leading colour. As demon- 



strated in 0, HI 



using the antenna subtraction scheme 
the explicit e-polcs in the dimension regularization pa- 
rameter of one- and two-loop matrix elements entering 
this calculation are cancelled in analytic and local form 
against the e-poles of the integrated antenna subtraction 
terms thereby enabling the computation of jet cross sec- 
tions at hadron colliders at NNLO accuracy. All of these 
techniques can be readily applied to the quark contribu- 
tions. 

For all of the observables considered here, we observed 
a dramatic reduction of the respective uncertainties in 
the theory prediction due to variations of the factoriza- 
tion and renormalization scales. We expect similar con- 
clusions when including the processes involving quarks. 
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